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FOREWORD 


This  is  the  final  report  furnished  by  Ultrusystems  {formerly 
Dynamic  Science)  on  the  "Analysis  of  a  Solid  Propellant  Gas  Generator 
using  an  Ammonium  Perchlorate  -  Polybutadiene  Composite  Propellant."  The 
report  was  prepared  for  the  Air  Force  Rocket  Pi  ''pulsion  Laboratory,  Edwards 
Air  Force  Base,  California.  It  is  submitted  in  accordance  with  the  provisions 
of  Air  Force  Contract  F0461I-71-C-0058  and  has  been  assigned  Dynamic 
Science  Report  Number  CS-72-5-1. 

Work  presented  in  this  report  began  in  June  1971  and  was  concluded 
in  Feb.  1972.  The  program  was  conducted  within  the  Engineering  and 
Computer  Sciences  Department  of  Dynamic  Science  at  their  Irvine  facility. 
Mr.  Douglas  Coats  was  principal  investigator  on  the  program  and  was 
assisted  by  Messrs.  G.  R.  Nickerson  and  H.  M.  Frey.  The  Air  Force 
Program  Monitor  was  Mr.  J.  Taska. 

This  technical  report  has  been  reviewed  and  is  approved. 


J.  Taska 

Project  Engineer,  AFRPL/LKCG 
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1.  INTRODUCTION 


The  purpose  of  this  contract  was  to  model  the  chemical  combustion 
and  mixing  process  in  the  Solid  Propellant  Gas  Generator  (SPGG)  device  and 
predict  experimental  conditions  wnich  were  most  favorable  for  laser  action 
in  the  exhaust  plume  of  the  device. 

The  basic  modeling  tool  to  be  used  was  Generalized  Kinetic  Analysis 
Program  (GKAP) ,  Reference  1,  and  the  tasks  involved  were  a)  the  establishment 
of  a  satisfactory  kinetic  combustion  mechanism  for  the  stream  tube  model, 

b)  determination  of  an  adequate  model  for  the  mixing  process  in  the  SPGG  device, 

c)  incorporation  of  the  chosen  models  into  the  computer  program,  and  d) 
comparison  between  analytical  predictions  and  experimental  observations 

The  propellant  system  in  the  SPGG  is  a  composite  made  up  of  crystaline 
Ammonium  Perchlorate  (AP)  in  a  Polybuiadiene  (PBAN)  fuel  matrix.  Much  of 
the  effort  under  the  contract  was  spent  determining  the  initial  composition  of 
AP  and  PBAN  in  the  gaseous  state  and  determining  a  reaction  set  which  would 
include  most  of  the  possible  mechanisms  for  the  combustion  process.  Section 
2  of  this  report  discusses  the  results  of  this  effort.  The  mixing  processes 
involved  in  propellants  of  this  type  were  also  extensively  investigated  and  the 
model  selected  for  the  SPGG  device  is  described  in  Section  3. 

In  order  to  compare  with  experimental  measurements  of  electromagnetic 
gain,  an  option  was  incorporated  into  the  GKAP  computer  code  to  calculate  tne  small 
signal  or  zero  power  gain  for  the  chemical  system.  The  method  of  calculation 
is  presented  in  Section  4. 

Section  5  reports  the  results  and  conclusions  arrived  at  from  the  total 
combustion  model  for  the  SPGG. 

During  the  course  of  the  contract,  experimental  data  obtained  by  the 
Air  Force  determined  that  the  SPGG  device  gave  negative  results  in  the 
measurement  of  gain.  For  this  reason  the  emphasis  on  the  work  performed 
under  this  contract  changed  from  a  detailed  analysis  of  the  AP/PBAN  propellant 
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system  of  the  SPGG  to  adding  a  more  complete  laser  modeling  capability 
to  the  GKAP  computer  program . 

Appendix  A  of  this  report  describes  two  mixing  models  incorporated 
into  the  GKAP  computer  code  during  this  contract  and  Appendix  B  shows  the 
result  of  a  H-F  ignition  study  to  verify  the  modeling  capability  of  GKAP, 
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2  .  CHEMICAL  MODEL  FOR  THE  AP-PBAN 
SOLID  PROPELLANT  SYSTEM 


The  AP-PBAN  solid  propellant  under  consideration  here  consists  of 
particles  of  crystalline  ammonium  perchlorate  (AP)  held  in  a  matrix  of  poly¬ 
butadiene  (PBAN)  binder.  The  weight  ratio  of  AP  to  PBAN  is  approximately 
70/30.  Unlike  rocket  engine  propellants,  no  aluminum  metal  is  present. 
During  the  combustion  process,  the  AP  crystals  are  thought  to  act  as  a 
monopropellant  producing  a  decomposition  flame,  which  in  turn  pyrolizes 
the  binder.  The  AP  decomposition  products  and  the  products  of  the  PBAN 
pyrolysis  then  react  exothermically  producing  further  combustion  products. 

A  flame  model  for  this  process  is  presented  in  Section  3. 

In  order  to  select  a  chemical  model  suitable  for  analyzing  the 
AP/PBAN  system,  the  combustion  process  has  therefore  been  considered 
to  occur  in  three  successive  steps.  These  steps  are: 

Step  1)  AP  decomposition 

Step  2)  PBAN  pyrolysis 

Step  3)  combustion  of  AP  decomposition  products 
with  the  PBAN  pyrolysis  products 

The  purpose  of  this  study  is  to  attempt  to  predict  by  use  of  a  kinetic 
model  which  gas  phase  species  will  result  from  the  combustion  process, 
and  by  what  reaction  mechanism  will  these  species  be  produced.  It  is  of 
particular  interest  to  be  able  to  identify  reactions  likely  to  take  place  which 
are  highly  exothermic  and  result  in  stable  species.  Species  formed  in  this 
way  will  tend  to  populate  excited  states  and  therefore  offer  the  possibility 
of  either  lasing  directly  or  of  transferring  their  energy  to  a  species  known 
to  lase  (such  as  CC^h  In  order  to  establish  a  reaction  mechanism  for  the 
AP/POAN  system  it  is  necessary  initially  to  consider  the  ground  state 
chemistry.  The  results  of  the  ground  state  analysis  can  then  be  used  to 
indicate  which  excited  state  species  are  probable. 


The  chemical  kinetics  of  combustion  for  a  composite  solid  propellant 
are  unknown  at  the  present  time.  Solid  phase  reactions  occurring  at  low 
temperature  below  the  burning  surface  may  be  of  importance  in  the  initial 
stages  of  combustion.  Reactions  of  this  type  are  known  to  be  compiex 
since  it  is  possible  for  many  kinds  of  relatively  stable  molecules  of  in¬ 
creased  molecular  weight  to  be  formed.  What  is  of  interest  in  this  study, 
are  the  gas  phase  reactions  v/hich  might  lead  to  excited  state  species. 
Consequently  a  simple  kinetic  model  has  been  adapted  for  the  initial  stages 
of  combustion,  which  considers  only  gas  phase  reactions.  The  subsequent 
stages  of  combustion  are  assumed  to  occur  at  temperatures  in  the  range 
known  to  exist  in  the  AP/PBAN  flame  (the  equilibrium  flame  temperature  at 
2.5  atm  is  2219°  K  for  the  AP/PBAN  used  in  this  study). 

The  combustion  chemistry  of  AP/PBAN  composite  propellant  burning 
is  only  partly  known  at  the  present  time.  Accurate  thermochemical  data 
(e.g.,  JANNAF  data)  is  available  for  most  chemical  species  of  interest,  but 
the  controlling  chemical  reactions  are  not  known.  It  is  possible,  however, 
to  postulate  a  system  of  chemical  reactions  sufficiently  comprehensive  to 
include  many  of  the  most  important  governing  steps  of  the  combustion  process. 
The  actual  steps  by  which  the  postulated  reactions  proceed  are  open  to 
question.  In  some  instances  excellent  experimental  reaction  rate  data  is 
available,  but  for  other  reactions  only  theoretical  estimates  are  available. 

Still  other  reactions  have  measured  data  available  v/hich  is  of  unknown 
accuracy. 

In  attempting  to  compile  a  reaction  mechanism  for  the  combustion 
process,  extensive  use  was  made  of  the  available  literature.  The  literature 
on  this  subject,  particularly  concerning  AP  combustion,  is  massive.  Never¬ 
theless  only  part  of  the  gas  phase  combustion  chemistry  for  this  system 
can  be  obtained  from  the  literature.  In  order  to  fill  in  the  missing  parts  of 
+he  reaction  chemistry,  reactions  have  been  postulated  which  appear  to  oe 
thermodynamically  plausible.  First  the  species  to  be  considered  were 
postulated.  These  are  presented  in  Table  2-1.  Certain  reactions  were  then 


2-2 


eliminated  due  to  steric  and  structural  considerations,  i.e.,  reactions 
involving  the  breaking  and  formation  of  several  or  bonds.  Some  reactions 
of  this  type  (i.e. ,  requiring  much  bond  rearrangement)  were  retained  because 
other  investigators  considered  them  important.  Bond  strengths  were  estimated 
using  the  JANNAF  data.  The  bond  dissociation  energy,  D,  for  molecule 

BC 


was  taken  as 


D  =  (AHf°  of  B)  +  (AHf°  of  C)  -  (aH£°  of  BC) 

where  the  heats  of  formation  are  at  298. 15°  K.  These  A  are  given  in 

Table  2-1  for  the  species  used  in  this  study.  The  species  are  listed  in 
alpha  numeric  order.  A  number  of  bond  dissociation  energies  estimated  in 
this  way  are  presented  in  Table  2-2. 

Particular  attention  was  given  so  as  to  include  reactions  occurring 
between  those  species  believed  to  be  prevalent  in  the  gas  phase  products. 
Reactions  which  occur  by  bond  breakage  through  collision  and  involve 
minimal  bond  rearrangement  were  given  high  priority,  regardless  of  whether 
the  reaction  was  endothermic  or  exothermic.  In  other  words  emphasis  was 
given  to  gas  phase  reactions  at  temperatures  near  the  flame  temperature. 

Where  it  was  necessary  to  estimate  reaction  rates,  order  of  magni¬ 
tude  estimates  were  used  .  There  are  of  course  many  methods 
by  which  reaction  rates  can  be  estimated.  For  example  bond-energy-bond 
methods,  collision  theory  methods,  the  Johnson  and  Parr  method  (Reference  2) 
and  the  Evans  and  Polanyi  method  (Reference  3  )  have  all  been  used  by 
various  researchers. 

The  methods  used  in  this  study  whenever  reasonably  reliable  rate 
data  was  not  available  are  presented  below. 


TABLE  2-2 


BOX’D  DISSOCIATION  ENERGIES 
ESTIMATES  FROM  JANNAF 
(D  kca  I/mole) 


dch-o 

= 

230.6 

= 

DCH3-H 

= 

101.9 

d:i-no 

= 

Dci-Cl 

57.3 

dh-o2 

Dci-o 

— 

64.3 

dx-h 

= 

Dcio-o 

= 

59 

dnh-h 

= 

■’co2-o 

= 

48 

dnh2-h 

= 

dcio-h 

- 

98.3 

dx-k 

= 

dci-oh 

= 

55.9 

a 

i 

o 

r 

XA-CO 

— 

54.7 

Dxo-o 

= 

5h-ci 

= 

103.1 

dX'-o2 

= 

5hc-ch 

rr 

2*29.8 

V° 

- 

\c-ch2 

' 

177.5 

do-h 

= 

*hci-h 

= 

111.9 

_ 

85.9 

91 

103.4 

226.0 

151 
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33.9 

40 

102.3 

119.2 


Order  of  Magnitude  Arrhenius  Equation  Rate  Estimates 


When  rate  constants  are  represented  by  an  Arrhenius  equation, 
k  =  AT  exp  (-E/RT) ,  where  T  is  the  absolute  temperature,  R  is  the  gas 
constant,  A  is  the  frequency  factor,  E  is  the  activation  energy,  and  n 
determines  the  pre -exponential  temperature  dependence,  order  of  magnitude 
approximations  have  been  made  as  follows: 

Exothermic  .  trimolecular  reactions 

B+C+M  ?  BC+M  ,  k  =  3xl016  t“°‘5 

Exothermic ,  bimoiecular  reactions  with  triatomic  transition  states. 

B+CD  t  BC+D  ,  k  =  5xlOU  T0,5  exp  -j-E/RT^ 

where  E  =  5.5%  of  the  CD  bond 
energy  (Hirschfelder  Rule, 
Reference  4  ) 

Exothermic ,  bimoiecular  reactions  with  transition  states  of  more 
than  three  atoms 


BC+DE  ?♦  BCD+E  ,  k  =  lxlOU  T°  *5  exp  -E/RT 

V.  j 

where  E  =  5.5%  of  the  DE  bond 
energy 

Exothermic  bimoiecular,  binary  exchange  reactions 

BC+DE  2  BD+CE  ,  k  =  lxlO10  T°  ’ 5  exp  {-E/RT 

where  E  =  28%  of  the  sum  of  the 
BC  and  DE  bond  energies 

The  reaction  rates  for  the  "spin  forbidden"  reactions  can  be  similarly  esti¬ 
mated  if  the  rate  constants  are  corrected  by  Boltzmann  factors  for  the  fact 
that  these  reactions  do  not  occur  in  the  ground  state . 

The  above  method  of  estimating  exothermic  reaction  rates  are  similar 
to  those  used  by  Tunder,  et  al,  Reference  5  . 
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AP  Decomposition 

>  ■ , 

The  Decomposition  anci  Combustion  of  Ammonium  Perchlorate- has  . 
been  the  subject  of  an  enormous  amount  of  research  during  the  past  20  years. 
A  systematic  review  of  this  subject  through  1968  is  given  in  Reference  6  . 
The  chemistry  of  this  process  must  still  be  regarded  as  unknown  to  a  large 
extent.  Under  conditions  of  interest  in  this  study  (steady  state  combustion 

X 

at  a  pressure  of  several  atmospheres)  it  has  been  determined,  however, 
that  an  AP  crystal  initially  undergoes  dissociative  sublimation  consistent 
with  the  reaction: 


NH*  C1C>4  (s)  =  NH3  (g)  +  H  CIO^  (g) 

i 

If  the  above  reaction  is  assumed  to  occur  near  the  AP  crystal  surface,  ihe 
chemistry  for  the  AP  system  reduces  to  an  oxidation  scheme  for  ammonia 
and  for  perchloric  acid.  Perchloric  acid  is  known  to  be  highly  unstable  and 
at  the  temperatures  which  exist  near  the  surface  of  the  propellant  will 
rapidly  decompose  by  the  following  steps: 


H  CIO „  -*  CIO-  +  OH 
4  3 

cio3  -*  CIO  +  02 

Several  chain  reaction  mechanisms  initiated  by  the  above  steps  have  been 
given  in  the  literature  (for  example  References  6  ,  7,8,  and  9  ).  For 
the  purposes  of  this  study  it  has  been  assumed  simply  that  NH^CIO^  initially 
breaks  down  through  the  absorption  of  energy  into  equal  moles  of  NH3,  OH, 
CIO,  and  02  according  to  the  mechanism 


nh4  cio4 
nh3,  IIC104 
OH.  C103 

mo,  o2 
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With  these  assumptions  a  kinetic  model  for  the  monopropellant 
combustion  of  AP  consists  of  ammonia  oxidation  in  the  presence  of  CIO. 

The  species  which  have  been  selected  as  being  of  possible  importance  are* 

Cl 

Cl  HO 
Cl  O 

ci2 

H 

H  Cl 
ho2 

h2 

h2o 

N 

NH 

nh2 

nh3 

NO 

no2 

N2 

n2o 

o 

OH 


*These  species  are  these  listed  in  Table  2-1  excluding  species  with  carbon. 
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It  is  necessary  to  consider  which  chemical  reactions  between  these 
species  are  likely  to  be  the  most  important.  For  the  ammonia  species  the 
the  following  reactions  were  selected  as  being  of  possible  importance.* 

NH3+C1  ?  NH2+H  Cl  '  .335  kcal/mole) 

NH3+C10  **  NH2+C1  HO  (  5. 150  kcal/mole) 

NH-.+H  NH  +H  {-  .762  kcal/mole) 

O  Z  Ca 

NH3+0  ^  NH2+OH  (  1.213  kcal/mole) 

NH3+OH  4  NH2  +  H20  (-15.890  kcal/mole) 

NH3+02  ^  NH2+H02  (  56.340  kcal/mole) 

The  NH3+OH  reaction  above  is  especially  important  because  of 
the  presence  of  large  amounts  of  both  species. 

Many  other  reactions  involving  NH3  are  of  course  possible  but  have 
not  been  included  for  various  reasons.  For  example  the  reaction 

NH3+°2  **  HN0+H20  (-  23  kcal/mole) 

has  not  been  included  because  it  is  stearically  complicated,  involving 
breaking  three  bonds,  and  the  activation  energy  is  so  high  that  it  can 
not  be  an  important  step. 

The  reduction  of  NH3  is  followed  by  the  reduction  of  NH3  and  of  NH. 
Twenty  three  reactions  have  been  included  to  account  for  this  process  and 
these  are  listed  as  reactions  74  through  96  in  Table  2-4  at  the  end  of 
this  section.  Table  2-4  contains  the  complete  reaction  table  used  in  the 


*The  heat  of  reaction  is  taken  as  the  sum  of  the  AH0,  of  the  products 

1 2  9  8 

(RHS)  minus  the  sum  of  AH°f  of  the  reactants  (LHS).  Thus,  if  the  sign 

*298 

is  positive  then  the  reaction  is  endothermic  left  to  right,  while  if  the  sign 
is  negative  then  the  reaction  is  exothermic  left  to  right. 


study  for  the  AP/PBAN  system.  The  format  and  contents  of  Table  2-4  are 
discussed  at  the  end  of  this  section. 

Other  reactions  occurring  between  the  species  selected  for  the  AP 
system  are  listed  as  reactions  44  through  64,  67  through  73  and  101  through 
119. 

The  above  set  of  74  reactions  include  all  of  the  21  reactions 
selected  by  McHale  for  AP  combustion  in  Reference  10  with  one  exception. 
In  the  present  study  the  reaction 

CLO+CLO  ♦♦  02+2Cl  {+  9 . 5  kcal/mole) 


was  used  rather  than 

CLO+CLO  £  02+CL2  (-48 . 4  kcal/mole) 


since  this  reaction  has  been  studied  and  found  to  proceed  by  the  former  path. 


PBAN  Pyrolysis 


Fundamental  studies  on  the  pyrolysis  of  polybutadiene,  such  as 
presented  in  Reference  11 ,  indicate  that  polvhutadiene  will  pyrolyze  into 
the  same  products  as  occur  in  chemical  equilibrium.  Results  of  equilibrium 
calculations  at  prescribed  temperature  and  pressure  are  presented  in 
Table  2-3.  A  total  of  33  species  were  considered  in  these  calculations. 

It  can  be  seen  that  the  major  equilibrium  products  at  the  lower  temperatures 
(600°  K)  are  ethylene  and  acetylene.  At  higher  temperatures  {2000°  K) 
ethylene  tends  to  disappear  but  of  course  acetylene  is  very  stable  at  high 
temperature.  It  has  been  estimated  that  the  surface  temperature  for 
AP/PBAN  composite  propellants  is  approximately  600°  C  to  700°  C  regardless 
of  pressure  (Reference  12).  At  this  temperature  as  can  be  seen  from 
Table  2-3  the  combined  mole  fraction  of  acetylene  and  ethylene  represents 
90%  of  the  equilibrium  product  gases. 

Experimental  studies  such  as  presented  in  Reference  11  show  that 
the  pyrolysis  of  PBAhl  results  in  formation  in  the  following  products: 

vinylcyclohexane 

butadiene 

ethylene 

acetylene 

In  other  words  the  polybutadiene  tends  to  pyrolyze  into  its  dimer, 
monomer,  and  fragments  of  the  monomer.  It  is  easy  to  see  how  this  might 
occur.  The  idealized  chain  structure  for  polybutadiene  is 

<-  CH2  -  CH  =  CH  -  CH2  -)n 


2-11 


TABLE  2-3 


EQUILIBRIUM  PRODUCTS  FOR  PBAN  PYROLYSIS 
ELEMENTAL  COMPOSITION  C?  u?  H1()  ^  O  m  N  j  j  PRESSURE  =  2.5  ATM 


Species 

600°K 

1000°  K 

Mole  Fractions 

1500°  K 

2000°  K 

ch4 

.000417 

.086639 

. 110060 

.005003 

CO 

.030559 

.028724 

.023087 

.019766 

c2H 

0 

0 

0 

.001276 

C2H2 

.390947 

.481691 

.598472 

.607228 

C2H4 

.556874 

.377117 

.066418 

.003111 

HCN 

.000833 

.011151 

.019324 

.020908 

H2 

0 

.000716 

.  176510 

.338903 

N2 

.003028 

.013886 

.006035 

.002979 

2 


The  weakest  links  in  the  chain,  are  links  1  and  2 . 


link  1  =41.9  kcal/mole 
link  2  =81.8  kcal/mole 

Breakage  at  link  1  results  in  the  monomer.  Further  breaks  might  then  occur 
at  link  2  such  that 


C4H6"C2H2+2CH2<C2H2+C2H4 
The  result  is  thus  the  formation  of  acetylene 

HC  s  CH 


and  ethylene 


H  H 


\ 

/ 

c 

=  c 

/ 

\ 

H  H 


both  of  which  are  thermodynamically  very  stable. 

Consequently,  for  the  purposes  of  this  study  acetylene  and  ethylene 
were  taken  as  the  primary  products  of  PBAN  pyrolysis.  The  combustion 
chemistry  is  thus  primarily  the  oxidation  chemistry  of  these  two  species. 

In  order  to  obtain  kinetic  mechanisms  for  the  oxidation  process,  reference 
was  made  to  the  existing  literature  on  the  subject. 

For  acetylene  the  reaction  mechanism  and  reaction  rates  given  by  Williams 
and  Smith  in  Reference  13  were  selected.  This  scheme  is  presented  in 
Table  2-4  as  reactions  2,23,24,30,31,32,33,34,35. 

For  ethylene  the  reaction  mechanism  and  reaction  rates  given  by  Sorenson 
et.  al.  in  Reference  14  were  selected.  This  kinetic  mechanism  is  a  subset 
of  their  more  comprehensive  kinetic  mechanism  for  ethane.  This  scheme 
is  presented  in  Table  2-4  as  reactions  7,  25,20,36,  and  37. 
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Combustion  of  AP  Decomposition  Products  with  PBAN  Pyrolysis  Products 

The  complete  AP/PBAN  reaction  set  used  in  this  study  is  presented 
in  Table  2-4.  In  addition  to  the  reactions  discussed  previously,  this 
table  contains  reactions  between  AP  and  PBAN  products.  Also  included 
are  third  body  reactions  between  major  species. 

Reactions  given  in  Table  2-4  include  all  of  those  reactions  evaluated 
by  the  University  of  Leeds  in  References  15  (  except  for  those  reactions  in¬ 
volving  the  H2C>2  species).  The  species  involved  in  these  12  reactions 
are:  CO,  CO^  H,  H02,  H2,  H20,  N,  NO,  K02,  N  ,  N20,  O,  OH,  O  . 

Heats  of  reaction  at  298.15°K  for  the  reactions  listed  in  Table  2-4 
are  presented  in  Table  2-5. 

The  format  of  Table  2-4  is  as  follows. 

Each  reaction  is  defined  by  one  line  of  information.  The  information 
given  on  each  line  read  left  to  right  is: 

Reaction;  reactants  =  products 

Rate  (left  to  right);  k  =  A  T  ^  e 

with  units  of  cc,  mole,  sec.,  °K,  kcal 

Reference 
Reaction  Number 

The  reactions  have  been  divided  into  two  types:  those  with  an  implied 
third  body  and  those  with  no  implied  third  body. 

Implied  Third  Body  Reactions 

These  reactions  assume  a  third  body  to  be  involved  in  the  collision 
process.  For  example  reaction  1  is 

co  +  o+  m  =  m  +  co2 

where  M  represents  any  of  the  other  species  present  in  the  system.  In  this 
study  all  third  bodies  were  assumed  to  have  the  same  efficiency.  These 
reactions  are  listed  in  alphanumeric  order  by  the  species  to  be  dissociated. 

Reactions  With  No  Implied  Third  Body 

These  reactions  are  listed  in  alphanumeric  order  by  the  principle 
species  to  be  reduced. 
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TABLE  2.4.  AB/l'BAM  KEAUTlON  SET  (continued) 
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3.0  MIXING  AND  COMBUSTION  MODEL  FOR  AN  AMMONIUM 
PERCHLORATE  AND  POLYBUTADIENE  COMPOSITE  SOLID  PROPELLANT 


For  most  rocket  motors,  both  solid  and  liquid  propellants,  it  is  assumed 
that  the  combustion  products  in  the  combustion  chamber  are  in  chemical 
equilibrium.  Mixing  is  usually  assumed  to  be  complete  or  is  handled  by 
computing  an  O/F  distribution  across  the  chamber.  The  above  assumptions 
are  usually  adequate  for  predicting  rocket  engine  performance.  However, 
the  emphasis  in  the  SPGG  device  is  on  the  possibility  of  the  development  of 
substantial  departures  from  chemical  equilibirum  in  the  combustion  chamber 
which  could  also  propagate  to  the  exit  plane  of  tne  nozzle.  The  two  physical 
processes  which  determine  the  chemical  composition  in  the  chamber  are  the 
mixing  rate  between  the  fuel  and  oxidizer  and  the  rate  of  chemical  combustion 
(kinetics).  Both  of  these  processes  are  coupled  together  and  depend  on  the 
physical  characteristics  of  the  device  itself. 

A  schematic  diagram  of  the  SPGG  device  is  shown  in  Figure  3-1.  The 
device  operates  at  a  chamber  pressure  of  2 . 5  atmospheres  with  a  contraction 
ratio  of  53.3:1.  The  propellants  are  an  ammonium  perchlorate  (AP)  oxidizer 
(73%  by  weight)  in  a  polybutadiene  (PBAN)  fuel  matrix  (22%  by  weight).  The 
AP  crystal  size  runs  from  200  microns  to  400  microns. 

If  the  assumption  of  one-dimensional  flow  in  chemical  equilibrium  is  made, 
the  velocity  of  the  combustion  product  gases  is  37.8  ft/sec  at  the  propellant 
surface  and  the  stay  time  in  the  nozzle  is  3.2  milliseconds.  From  the  velocity 
and  equilibrium  composition  the  Reynolds  number  in  the  chamber  was  estimated 
to  be  on  the  order  of  5  to  10  based  on  AP  particle  sizes  of  200  to  400  microns, 
respectively.  Reynolds  numbers  in  this  range  preclude  turbulent  mixing  on  the 
scale  of  the  AP  particle  diameter  due  to  viscous  forces.  Photographic  data  taken 
of  the  AP-PBAN  propellant  combustion  process  in  the  SPGG  device  also  confirms 
the  lack  of  turbulence  on  this  scale.  The  motion  picture  show  well  established 
streamers  many  AP  particle  diameters  above  the  propellant  surface  and  are  similar  in 
nature  to  columnar  diffusion  flames.  However,  it  is  not  possible  to  determine  if  the 
visible  streamers  are  part  of  primary  or  secondary  combustion  zone  without  accurate 
temperature  measurements.  With  the  presence  of  apparent  columnar  diffusion 
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iqure  3-1.  SPGG  Device  Schematic 


flames,  it  is  tempting  to  assume  a  laminar  diffusion  flame  model  to  describe 
the  mixing  and  kinetics.  Unfortunately,  however,  there  are  two  main  draw¬ 
backs  to  this  approach,  the  first  is  that  the  columnar  diffusion  flame  theory 
predicts  that  the  height  of  the  flame  zone,  i.e.  ,  the  height  at  which  the  tempera¬ 
ture  of  the  gas  becomes  approximately  the  equilibrium  flame  temperature,  is 
independent  of  pressure,  which  is  in  contradiction  with  experimental  data  and 
conclusions  of  References  16,  17,  18,  19.  The  second  drawback  is  that 
there  is  no  a  priori  method  of  establishing  the  composition  and  energy 
content  of  the  fuel  and  oxidizer  streams.  It  becomes  apparent  then,  that 
any  rigorous  treatment  of  the  problem  of  composite  propellant  burning  must 
consider  the  interaction  between  the  combustion  flames  and  the  propellant 
surface. 

This  very  problem,  the  steady  state  combustion  of  composite  propellants, 
has  been  studied  extensively  for  years.  However,  despite  the  large  amounts 
of  experimental  and  theoretical  research  done  on  this  problem,  no  definitive 
model  has  been  advanced  which  is  capable  of  accurately  predicting  the 
combustion  characteristics  and  flame  structure  of  a  composite  solid  propellant. 

Most  of  the  models  proposed  to  date  Ref.  16,  17,  and  18  are  multi- 
flame  models  and  are  for  small  oxidizer  grains  (3  to  20  microns), 
moderate  to  high  pressures  (100  psia  and  greater) ,  and  are  intended  to  pre¬ 
dict  the  linear  burning  rate  of  the  composite  propellant.  However,  in  the 
SPGG  device,  we  have  noted  that  the  oxidizer  grain  size  is  quite  large  (200 
to  400  microns),  and  the  pressure  low  (2.5  atm),  so  that  the  applicability 
of  the  models  proposed  becomes  unclear.  Another  disadvantage  of  the  models 
proposed  in  the  literature  is  that  they  are  essentially  one-dimensional  in 
nature  and  ignore  or  gloss  over  the  details  of  the  mixing  process  between  the 
gaseous  fuel  and  oxidizer.  The  mixing  process,  however,  is  precisely  the 
point  of  interest  for  the  SPGG  device  since  we  are  looking  for  nonequilibrium 
effects  which  could  be  caused  by  mixing. 

In  addition  to  the  mixing,  the  source  of  energy  needed  to  perform  the 
endothermic  pyrolysis  of  the  fuel  binder  must  be  known  in  order  to  establish 
the  energy  content  of  the  fuel  and  oxidizer  in  the  gaseous  state.  The 
possible  sources  of  this  energy  by  way  of  conductive  heat  transfer  and 
combustion  are  shown  in  Figure  3-2  and  are: 
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Binder  Pyrolysis 


Figure  3-2.  Illustrative  Schematic  of  Combustion  Processes 
In  An  AP/PBAN  Composite  Solid  Propellant 
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1.  Exothermic  decomposition  of  the  oxidizer  crystal. 

2.  Exothermic  gas  phase  reactions  between  the  oxidizer 
decomposition  products  and  the  fuel  pyrolysis  products. 

3.  Exothermic  heterogeneous  reaction  between  the  binder 
and  the  oxidizer  interface. 

The  rate  controlling  step  for  the  combustion  process,  while  not  known 
exactly,  must  consist  of  one  or  more  of  the  following: 

a.  Chemical  reaction  kinetics  occurring  in  the  various 
reaction  zones . 

b.  Mixing  (diffusive  and/or  convective)  between  the 
oxidizer  decomposition  products  and  fuel  pyrolysis 
products . 

c.  Conductive  heat  transfer  from  the  flame  into  the  propellant 
surface . 

Unfortunately,  there  is  a  lack  of  agreement  (e.g.  Reference  19  and  20) 
in  the  results  of  experiments  designed  to  determine  the  structure  of  the 
flame  zone.  It  is  particularly  important  to  determine  by  measurement 
the  temperature  profile  above  and  below  the  propellant  surface.  This 
knowledge  is  necessary  in  order  to  determine  the  importance  of  solid 
phase  reactions  in  sustaining  the  combustion  process  and  to  answer 
questions  regarding  the  transfer  of  energy  from  gas  phase  chemical  re¬ 
actions  to  the  prepellant  surface.  Experimental  results  to  date  are  not 
sufficiently  accurate  to  provide  this  type  of  temperature  profile  data. 
These  results  do  show,  however,  that  temperatures  near  the  equilibrium 
flame  temperature  for  the  mixture  occur  very  close  to  the  propellant  sur¬ 
face. 


Reference  17  has  proposed  a  three  flame  model  for  composite  propellant 
burning.  The  model  consists  of  an  AP  deflagration  flame  extending  from  the 
surface,  a  primary  flame  between  the  AP  decomposition  products  and  the 
PBAN  pyrolysis  products,  and  a  final  or  secondary  diffusion  flame  between  the 
basically  oxidizer  and  fuel  streams.  This  model  would  fit  nicely  with  the 
observations  of  the  combustion  process  in  the  SPGG  device.  While  the  model 
is  essentially  one-dirnensional  with  the  two-dimensional  effects  accounted 
for  by  varying  the  time  constants  for  the  various  physical  mechanisms 
involved,  it  has  great  appeal  if  a  reasonable  method  could  be  established 
to  determine  the  enthalpy  and  composition  of  the  first  two  flames.  To 
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determine  if  a  reasonable  method  is  obtainable,  it  is  necessary  to  get  an 
estimate  on  the  amounts  of  energy  transferred. 

A  simple  equilibrium  model  was  used  to  determine  the  amount  of 
energy  required  to  pyrolyze  the  fuel  binder.  In  this  model  it  was  assumed 
that  the  gaseous  fuel  products  were  in  equilibrium  at  the  surface  temperature 
(also  pyrolysis  temperature)  of  the  fuel .  The  surface  temperature  of  the  fuel 
was  taken  as  600°C<  Pef.  11  ,  and  it  was  found  that  approximately  1250 
cal/gm  of  heat  addition  was  necessary  to  pyrolyze  the  fuel.  This  result 
compares  favorably  with  estimates  made  from  the  activation  energy  of  PBAN. 

The  same  type  of  calculation  was  performed  for  AP  and  it  was  found  that 
the  AP  equilibrium  products  at  the  surface  temperature  of  propellant  had 
an  enthalpy  of  173  cal/gm  lower  than  the  enthalpy  of  AP  solid.  This 
difference  in  enthalpy  can  be  explained  by  assuming  that  the  energy  difference 
would  be  used  to  pyrolyze  the  PBAN.  However,  using  the  relative  weight 
fractions  of  the  AP  and  PBAN  in  the  composite  propellant,  it  was  calculated 
that  the  AP  would  have  to  give  up  352.3  cal/gm  in  order  to  pyrolize 
the  PBAN.  Hence  it  becomes  evident  that  in  order  to  establish  the  energy 
content  of  the  two  streams  it  is  necessary  to  know  the  amount,  of  energy 
transferred  by  the  three  mechanisms  stated  before  and  where  in  the  gene¬ 
ration  of  the  two  streams  it  is  introduced  since  the  chemical  kinetics  are 
highly  sensitive  to  temperature  variations.  In  short,  to  model  accurately 
the  propellant  burning,  not  even  considering  solid  phase  reactions,  it 
would  be  necessary  to  do  a  complete  two-dimensional  calculation  with 
kinetics  and  heat  and  mass  transfer.  This  task  is  clearly  not  within 
the  scope  of  the  present  contract. 

An  alternative  to  a  complete  calculation  is  to  assume  that  the  overall 
combustion  process  is  kinetically  controlled  since  the  mixing  is  very  quick 
and  the  kinetics  are  slow  at  these  low  pressures.  This  assumption  is  standard 
for  low  pressure  systems  such  as  SPGG;  however,  the  very  large  grain  sizes 
of  the  AP  crystals  cast  some  doubt  on  the  validity  of  the  quick  mixing  part  of 
the  assumption.  If  this  assumption  is  accepted,  however,  then  we  can  assume 
that  the  decomposition  products  of  the  oxidizer  and  binder  as  discussed  in  Section 
2  are  injected  into  an  environment  of  a  premixed  flame  and  that  the  heat  transfer 
necessary  to  drive  the  combustion  is  instantaneous  so  that  the  enthalpy  of  the 
flame  is  that  of  the  composite  propellant.  The  composition  of  the  premixed  flame 
is  the  only  unknown  quantity  left  in  the  process.  We  shall  thus  make  the  further 
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assumption  that  the  premixed  flame  can  be  characterized  by  an  equivalent 
flame  in  a  well  stirred  reactor  whose  length  and  stay  time  are  the  same  as 
those  in  the  overall  mixing  region  in  the  combustion  chamber.  The  only 
justification  for  the  last  assumption  is  that  it  will  allow  a  hopefully  reason¬ 
able  nonequlibrium  chemical  result  with  energy  being  conserved. 

In  order  to  perform  the  calculation  as  stated  the  overall  mixing  length 
which  includes  both  the  primary  and  secondary  diffusion  flames  must  be 
known.  To  a  first  order  approximation  the  overall  mixing  length  is 

*  ~  Tvs  (3-D 

where  r  =  characteristic  time  for  mixing 

v  =  velocity  normal  to  the  surface 
of  the  propellant 

r  is  defined  for  molecular  diffusion  and  convection  as: 

T  =  (diffusion)  (3-2) 

r  =  7T-  (convection)  (3 -2b) 

VT 

where  d  =  distance  to  mix  through 

£  =  binary  diffusion  coefficients 

Vj  =  velocity  tangent  to  the  propellant 
surface 

Making  use  of  kinetic  theory  and  noting  that  d  =  1/2 ft,  where  6  is  equal  to 
the  particle  diameter,  we  get: 

62 (P/P  )v 

^  =  - °i  S7£  (diffusion)  (3-3a) 

■^(T/Tor- /fa 


(convection) 


(3  —3b) 


The  mixing  length  for  pure  diffusion  is  easily  calculated  to  be  24.6  6  for 
400  micron  particles  and  6.  15  6  for  2  00  micron  particles.  However,  in  order 
to  calculate  the  convective  mixing  length  we  must  estimate  the  tagential  velocity, 
Vrp.  Due  to  the  surface  roughness  and  irregular  erosion  rate  the  tangential 
velocity  should  be  on  the  order  of  5  to  10%  of  the  one-dimensional  normal 
velocity.  This  estimate  would  yield  a  mixing  length  of  56  to  106.  It  is 
obvious  then  that  the  diffusional  mixing  length  will  be  the  maximum  possible 
for  the  combustor. 

Equation  (3-1)  is  of  course  only  a  rough  approximation  to  the  real 
mixing  length  since  it  completely  ignores  variations  in  concentration  and 
temperature  gradients  caused  by  the  mixing  and  combustion.  An  exact 
solution  of  the  multi-component  diffusion  equation  is  not  within  the  scope 
of  this  project,  however,  we  can  consider  the  simpler  case  of  binary 
diffusion.  From  Fick's  Law  the  binary  diffusion  equation  without  convection 
is: 


S  x: 

p  ^  '  (p£  7x.)  (3-4) 

where  o  =  density 

th 

x.  =  mass  fraction  of  the  r  species 
Jfr  =  binary  diffusion  coefficient 


Noting  the  similarity  between  mas:'  n  if  fusion  ar.ci  :nr- rr:.ui  diffusion, 

equation  (3-4)  can  be  solved  numerically  using  a  lumDea  oarameter  electrical 
analog  finite  difference  method  as  is  done  in  heat  transfer  problems,  see 
Ref.  21  .  This  method  assumes  that  the  problem  can  be  treated  as  a  finite 
network  of  nodes  with  each  node  having  uniform  properties.  Under  these 
assumptions  equation  (3-4)  reduces  to 


x. 

_L 


-  x. 

i 


) 


(3-5) 


.  1_ 

where  Cl  =  capacitance  of  the  1  node 

R.j  =  resistance  between  node  i  and  j 


It  can  be  shown  tnat 


Rij  pAJt 
C.  =  pAir 

cross-sectional  area 
distance  between  node  i  and  j 

length  of  node  i 


where  A 


Ar..  = 


ari  = 


(3-6) 


If  we  only  consider  a  simple  2  node  system  with  constant  properties,  then 
equation  (3-5)  reduces  to 


dx. 


12  1 


((l-x,)o  -x,) 


(3-7) 


where  0  =  voIume  of  node  1 
volume  of  node  2 


The  solution  of  (3-7)  with  the  boundary  conditions 


x1(t=0)  =  1 
x2(t  =  0)  =  0 


is 


x,  = 


a  +  e 
l  +o 


-  (l  +o)t 
RC 


(3-8) 
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The  RC  product  for  the  2  node  system  becomes 


RC  = 


£r 


12’  &r 
3 


1 


Figure  3-3  shows  the  results  of  a  computation  made  using  the  two 
node  model  for  a  400 n  diameter  particle  assuming  a  diffusion  coefficient 
of  1 .  As  can  be  seen  from  the  figure,  there  is  an  appreciable  amount  of 
mixing  within  the  first  particle  diameter  and  mixing  is  complete  soon  after 
6  particle  diameters. 

The  two  node  model  predicts  a  higher  initial  rate  of  mixing  than  will 
occur  in  reality  due  to  the  high  initial  concentration  gradients.  A  multi-node 
model  computation  was  attempted  by  using  a  computer  program  to  solve 
equation  (3-5)  under  the  following  boundary  conditions 


S  x . 
_ i 

S  r 


-  0 

T  =  0 


r  rfuel 


=  0 


and  assumptions 


_S_ 

S  r 


>>> 


_8_ 
a  e 


where  r  =  radial  direction 
z  =  normal  direction 
9  =  azimuthal  direction 


with  density  being  calculated  by  assuming  constant  molecular  concentration 
at  each  node  and  temperature  evaluated  from  an  equilibrium  table  of  O/F 
versus  temperature.  The  results  of  the  computation  were  not  entirely 
satisfactory  since  a  very  high  mixing  rate  was  initially  predicted  at  the 
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OX  IN  OXIDIZER  STREAM 


OX  IN  FUEL  STREAM 


PARTICLE  DIAMETERS 


Figure  3-2.  RESULTS  OF  2  NODE  DIFFUSION  MIXING  MODEL 
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juncture  os  the  two  streets  which  negates  the  assuctptiaa  that  the  deriwatires 
:..  the  z  direction  <*re  ssill.  Howcw.  the  model  did  ^cdc;  that  sat 
auxin;  is  completed  an  »s;e  of  ZtKOjs  to  depeodiag  a®  the  A?  picticSt 

size  used. 

The  above  results  for  the  mixing  length  were  used  so  mifce  two 
stirred  reactor  runs  with  .^ettar  lengths  at  2409®  a~d  SSSO®  w£ath  w ere 
the  licittig  lengths  calculated.  The  results  of  these  stored  reactar  «kd*iau 
were  then  used  as  inputs  into  a  streantoce  caJcsslata-cca  for  the  chamiiw 
and  nozzle  of  the  S?GG  device.  The  results  off  these  conrguta  tines  arc- 
presented  in  Section  5  of  this  report. 

Since  the  stirred  reactor  model  was  used  to  predict  the  nooe^uatiirim 
start  conditions  in  the  chamber  off  the  SPGG  device,  a  brief  desosptjoei  off 
the  node!  is  presented  in  Appendix  A.  Also  presented  cs  Appendix  A  is  a 
description  of  a  shear  layer  mixing  option  added  to  the  -C&A?  co erp3S.tr 
prograa  under  this  contract. 

In  conclusion,  the  fixing  and  combustion  off  the  composite  pcope.lana 
used  in  the  SPGG  device  was  studied  extensively.  The  literature  was  searched 
and  it  was  found  that  no  satisfactory  model  of  this  phenomena  existed 
for  the  purpose  of  determining  propellant  flame  structure  and  cocposilioc- 
Various  models  were  tried,  including  a  pure  diffusion  model,  with  little 
success  due  to  the  difficulty  and  two  dimensional  nature  of  the  physical 
problem.  In  the  end  a  stirred  reactor  model  was  adopted  since  this  mode! 
allowed  for  finite  rate  chemistry  'while  conserving  energy  between  'he 
fuel  and  oxidizer. 


4.  SMALL  SIGNAL  GAIN  CALCULATION 


The  primary  purpose  of  the  SPGG  device  was  to  determine  if  laseable 
species  could  be  produced  in  quantity  in  the  exhaust  plume  and  the  purpose 
of  this  contract  was  to  analyze  the  combustion  process  in  the  device  and 
determine  the  mechanism,  if  any,  which  would  produce  the  laseable  species. 


Experimentally  the  presence  of  laseable  species  is  detected  by 
measuring  the  gain  or  loss  in  intensity  of  a  C02  laser  beam  which  penetrates 
the  vibrationally  excited  media.  In  order  to  correlate  with  experimental  data 
and  also  since  the  presence  of  vibrationally  excited  species  does  not 
insure  laseability  it  was  necessary  to  include  a  small  signal  gain  calculation 
option  into  the  GKA.P  computer  code. 

The  definition  of  the  small  gain  coefficient,  Qtg,  is  that  if  a  plane 
wove  of  intensity  is  propagating  in  the  positive  x  direction,  at  x  =  Xq, 
then  the  increase  or  decrease  in  intensity  with  x  is  given  by 


I 


(x-x0) 


(4-1) 


the  subscript  0  is  appended  to  a  to  denote  this  result  as  the  "small  signal 
gain,"  i .  e .  ,  the  limiting  value  for  I  -*  0  . 

The  equation  for  the  gain  of  electromagnetic  radiation  was  originally 
derived  (independently)  by  Fuchtbauer  and  Fadenberg  and  was  intended  to 
describe  absorption.  From  Refs.  22,  23,  and  24  the  gain  coefficient  is  given 

by 


a. 


A  c2 
mn 

8  rrp2 


(4-2) 


The  amplification  is  associated  with  transition  between  two  vibrational  states,  v' 
and  v.  The  prime  denotes  the  upper  energy  level.  The  radiation  emitted 
by  the  vibrational  transfer  in  energy  levels  is  not  monochromatic,  but  is 
spread  out  in  frequency  from  a  combination  of  Doppler  broadening,  due  to 
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:=*o2ccular  solsons.  Lc?rer.:z  broadening,  due  to  collisions,  and  natural 

broadening,  due  to  radiation  decay.  The  frequency  distribution  may  be 

represented  by  a  line  shape  function,  G  h/.v  ),  v/hich  has  roughly  a  bell 

ran 

shape.  Of  ~ost  interest  is  the  line  function  at  the  line  center  since  it  is 
a  rsaadma  there.  For  natural  and  collisional  (Lorenz)  broadening  at  the  line 
center. 


G,J'-“/r,r.)  =  77 1 


(4-3) 


where  (;  v)  ,  /,  is  the  full  line  width  at  half  maximum  and  is  equal  to 


v)  for  natural  broadening 

(y  =  radiation  dampening  factor) 


i  2  =  —  for  Lorenz  broadening 

(Z  =  collision  frequency) 


•"or  Doppler  broadening  the  line  function  at  the  line  center  is 


G<"’I'm„> 


/ 1  n2  _ 2_ 

V  7T  (WD) 


(4-4) 


where  («LV_J  is  given  by 


(LV  _J  =  2  '1q-2  v  Vt/M 


Not  included  in  equation  (4-2)  or  considered  here  is  the  variation  of 
the  frequency  with  rotational  energy  levels,  J  quantum  numbers,  and  the 
vibrational-rotational  interaction  factor.  These  effects  are  of  second  order. 


The  Einstein  A  coefficient,  A  ,  is  related  to  the  inverse  of  the  decay  time, 

mn 

r,  for  the  transition  from  the  upper  energy  level,  (v‘ ,  Jm)  to  the  lower  energy 
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level,  (v,  Jn) ,  Since  the  only  J  dependent  part  of  Amn  considered  here  is 

the  dipole  moment  of  the  molecule,  then  A  may  be  written  as  the  product 

mn 

of  J  dependent  and  independent  parts  by  introducing  the  line  strength,  Sj  , 
that  is,  n 


A 


mn 


(4-5) 


where  S.  =  J  for  P-branch  transitions 

Jn 

”  Jn+1  for  R-branch  transitions 
v1  1 

Av  =  — —  =  Einstein  probability  of 
Tv  spontaneous  emission 

corresponding  to  the 
integrated  band  intensity 
(sec-  M 


Hence  equation  (4-2)  may  be  rewritten  as 


% 


Av  A  : 
v 

8  77 


N  , 

V 


m 


N 


3m 


g 

3  i  — * 


mn 


(4-6) 


The  degeneracies  or  statistical  weights,  gm  and  gn,  are  given  by 

gk  =  (2Jk  +  D  (4-7) 

since  the  freedom  to  rotate  about  any  axis  allows  2  J  +  l  states  to  have  the 
same  energy  for  a  given  value  of  J. 

The  total  number  densities  of  the  vibrationally  excited  species,  N  ‘s, 
are  calculated  directly  by  the  GKAP  computer  program  using  finite  rate 
chemistry.  However,  the  gain  equation  requires  that  the  distribution  over 
the  rotational  energy  levels  for  each  vibrational  level  be  known.  Since  it 
is  well-known  that  the  rotational  degrees  of  freedom  in  a  gas  equilibrate 
rapidly  among  themselves  and  with  the  translational  degrees  of  freedom, 
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one  can  reasonably  assume  that  the  molecules  are  distributed  over  the 
rotational  states  according  to  the  Boltzmann  distribution  with  rotational 
temperature,  T  ,  equal  to  the  translational  temperature,  i.e,,  the  "gas 
temperature."  With  this  assumption  the  number  densities,  N  T ,  within  a 
vibrational  state  for  a  given  J  are 


N 


v,J 


N  B  he 

q-2-  9j  EXP  (-1(1  +  1)  ) 


where  N 

v 


total  number  density  of  the 
vibrational  state  (molecules/cc) 

rotational  constant  (cm--*-) 

Planck's  constant  (erg-sec) 
Boltzmann's  constant  (erg/°K) 
speed  of  light  (cm/sec) 
rotational  partition  function 


(4-8) 


The  rotational  partition  function  for  small  B  and  large  T  has  been  shown 
to  be 

^ro,  *  <^>/6 

where  6  =  symmetry  factor  (see  Ref.  25  and  26) 


For  linear  symmetric  triatomic  molecules  such  as  CO^  the  symmetry 
factor  is  2  and  reflects  the  selection  rule  for  such  molecules.  The  selection 
rule  being  that  for  the  rotational  states  only  even  or  odd  values  of  J  are 
allowed  for  a  given  vibrational  state.  The  net  effect  is  that  the  populations 
for  a  given  J  are  doubled  and  hence  the  small  signal  gain  coefficient  is 
also  doubled . 

As  implemented  in  the  computer  program  the  variation  of  B  with  v 
in  equation  (4-9)  is  ignored  since  the  variation  is  small. 

Thus,  equations  (4-3)  or  (4-4),  (4-6),  (4-7),  (4-8),  and  (4-9)  are 
sufficient  for  the  computer  program  to  calculate  the  gain  coefficient  given 
the  type  of  broadening ,  A^,  X,  B,  and  6  for  a  given  J.  The  program 
calculates  and  prints  out  the  maximum  value  of  for  hoth  P-branch  and 
R-branch  transitions  for  a  given  range  of  J's. 


5.  RESULTS  AND  CONCLUSIONS 


As  stated  previously,  the  purpose  of  this  effort  was  to  analyze  the 
combustion  process  in  the  SPGG  device  to  find  out  the  mechanism,  if  any, 
which  could  produce  quantities  of  excited  state  species  sufficient  to  lase. 
However,  experimental  results  obtained  by  the  Air  Force  during  the  course 
of  the  contract  showed  completely  negative  results  in  the  measurements  of 
gain.  For  this  reason  the  emphasis  on  the  work  performed  under  this  contract 
changed  from  a  detailed  analysis  of  the  AP/PBAN  propellant  system  of  the 
SPGG  device  to  adding  a  more  general  modeling  and  analysis  capability  to 
the  GKAP  computer  code.  Thus  the  analysis  of  the  SPGG  device  is  not  totally 
complete,  however,  it  does  show  some  interesting  results  and  trends. 

The  mixing  and  combustion  model  selected  for  tne  SPGG  device, 
see  Section  3,  was  that  of  a  stirred  reactor  of  a  length  of  the  overall  mixing 
zone.  The  input  to  the  reactor  consisted  of  the  decomposition  products  of 
the  AP  and  the  pyrolysis  products  of  the  PBAN  as  discussed  in  Section  2.  The 
chemical  kinetic  reaction  scheme  of  that  section  was  also  used  as  the  mecha¬ 
nism  within  the  reactor. 

Two  stirred  reactor  runs  were  made  using  mixing  lengths  of  2400/t 
(6  AP  particle  diameters)  and  9600/2  (24  diameters).  The  results  of  these 
calculations  and  a  chemical  equilibrium  calculation  are  shown  in  Table  5-1 
for  the  major  thermodynamically  stable  species. 


TABLE  5-1 

SPECIE  CONCENTRATION 
FOR  COMBUSTION  MODELS 


Calculation 

SPECIE  (MOLE  FRACTIONS) 

CO 

O 

O 

H  Cl 

H2 

h2o 

N2 

Chemical 

Equilibrium 

.2868 

.04006 

. 138321 

.264402 

.194377 

.073011 

Stirred 

Reactor 
(2400  u) 

.1897 

.0727 

. 150575 

.1567 

.2770 

.0748 

Stirred 

Reactor 
(9600  p) 

.20748 

.06518 

.  1497 

.  1704 

.2707 

.07507 
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This  table  shows  a  substantial  departure  from  equilibrium  for  all 
species  for  both  the  6  and  24  diameter  cases.  As  is  expected  from  a  stirred 
reactor  model,  the  longer  the  length  (stay  time),  the  closer  to  equilibrium 
the  output  composition  will  be. 

Examining  the  reasons  for  the  non-equilibrium  composition,  the 
following  becomes  clear  for  the  species  listed  below. 

C02  and  CO 

The  reaction  CO+CIO-^CO^+CI  is  the  dominant  mechanism  for  the 
overproduction  of  C02  and  deficiency  of  CO.  The  large  excess  of  CIO  being 
introduced  to  the  reactor  causes  this  reaction  to  proceed  much  faster  than 
the  other  mechanism  provided  for  destroying  C02  and  producing  CO. 

HC1 

The  overproduction  of  HC1  is  caused  by  the  reaction  C1+H2-»HC1+H 
which  is  producing  HC1  at  a  rate  3  times  faster  than  the  largest  HC1  oxidation 
reaction,  HC1+0-+C1+0H ,  is  destroying  it.  The  Cl  for  this  reaction  is  being 
produced  by  C0+C10-+C02+C1  reaction  and  the  by  the  NH^  reduction  reactions, 
NH3+H-*NH2+H2  and  NH2+H-»NH+H2 .  Note  that  CO,  CIO,  and  NH3  are  input 
reactants  to  the  reactor. 


The  underproduction  of  H2  is  due  mainly  to  the  overproduction  of  HC1 
since  most  of  the  H2  produced  by  the  reduction  of  NH^  is  consumed  by  the 
production  of  HC1. 

H2° 

The  excess  amounts  of  H20  are  produced  by  the  following  reactions 
with  the  OH  radical  which  is  an  input  reactant  to  the  reactor 

ch2o+oh  -4  h2o+hco 
h2+oh  -»  h2o+h 

NH+OH  -4  H20+N 
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The  formaldehyde  for  the  first  reaction  comes  from  the  oxidation  of  acetylene 


c2h2+oh  •♦ch2o+ch3 


and  the  H2  and  NH  are  produced  by  the  reduction  of  ammonia.  Both  C2H2 
and  NH^  are  present  in  large  quantities  since  they  are  inputs  to  the  reactor. 


The  slight  excess  amounts  of  N2  above  the  equilibrium  concentration 
is  due  the  fact  that  N2  is  added  directly  into  the  reactor  and  that  there  are  no 
significant  reactions  destroying  N2  at  these  temperatures. 

Thus  from  the  above  analysis  the  departure  from  equilibrium  can  be 
attributed  to  the  presence  of  large  quantities  of  the  raw  reactants  which  do 
not  have  sufficient  time  to  react  at  the  chamber  pressure  of  the  SPGG  device. 
For  summary  purposes  the  input  reactants  and  the  major  reaction  mechanisms 
within  the  reactor  are  shown  below. 

REACTOR  KINETIC  MECHANISM 


INPUT  REACTANTS 


OUTPUT  CONCENTRATIONS 


NH3 

CO+CIO 

C02+C1 

CIO 

Cl+H2 

HC1+H 

°2 

NH3+H 

-4 

nh2+h2 

OH 

nh2+h 

nh+h2 

C?H? 

ch2o+oh 

h2o+hco 

o 

PO 

h2+oh 

H20+H 

CO 

NH+OH 

h2o+n 

N2 

c2h2+oh 

-4 

ch2ok;h3 

Greater  than  equilibrium 

co2 

HC1 

h20 


Less  than  equilibrium 
CO 


Figure  5-1 .  STIRRED  REACTOR  CHEMISTRY 
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Kinetic  stream  tube  runs  were  made  for  the  nozzle  with  a  variety  of 
initial  conditions  ranging  from  full  chemical  equilibrium  start  conditions  to 
the  output  of  the  stirred  reactor.  Departures  from  chemical  equilibrium  at 
the  throat  plane  of  the  nozzle  were  noted  in  all  cases  and  varied  from  small 
to  large. 

In  those  runs  started  from  equilibrium  or  close  to  equilibrium,  the 
flow  stayed  in  or  equilibrated  rapidly  until  the  vicinity  of  the  throat  plane. 

At  approximately  an  area  ratio  of  2  the  flow  experiences  a  very  rapid  accele¬ 
ration  from  a  Mach  number  of  M=.25  to  M=1.0  in  the  distance  of  0.5  inches. 
The  high  velocities  and  corresponding  short  stay  times  in  this  region  does 
not  allow  the  chemistry  to  remain  in  equilibrium  at  the  low  pressures  as  are 
characteristic  in  the  SPCG  device. 

Cases  which  were  run  from  stirred  reactor  initial  conditions  stayed 
substantially  out  of  equilibrium  for  the  entire  time  in  the  nozzle.  Table  5-2 
lists  selected  specie  concentrations  at  the  throat  (also  exit)  plane  of  the 
nozzle  for  the  various  cases  of  interest. 


TABLE  5-2 

SPECIE  CONCENTRATIONS  AT  THE 
NOZZLE  EXIT  PLANE 


CASE 

SPECIE  CONCENTRATIONS  (MOLE  FRACTIONS) 

Full  Chem. 

Eq. 

CO 

CM 

o 

o 

HC1 

H2 

h20 

N 

2 

.2824 

.04483 

.139018 

.26976 

.  190044 

.073073 

Chem.  Eq. 

Start 

.2864 

.0404 

.1383 

.2647 

.194096 

.073015 

Stirred  Reactor 
Start 

(L=2400m) 

.2300 

.05455 

.14664 

.20136 

.2549 

.07707 

The  results  of  these  calculations  show  that  any  significant  departure 
from  chemical  equilibrium  will  persist  throughout  the  chamber  of  the  device. 

The  most  likely  reason  for  this  effect  is  the  low  chamber  pressure  of  the 
device  since  the  reaction  rates  proceed  approximately  as  Pff  where  a  is  the 
order  of  the  reaction . 

However,  departure  from  chemical  equilibrium  does  not  necessarily 
mean  that  a  population  inversion  has  occurred  among  the  vibrationally  excited 
states  of  a  species .  For  this  reason  gain  calculations  were  performed  on  C02 
for  the  L=2400/x  stirred  reactor  start  case.  The  excited  state  C02  reaction 
set  is  shown  in  Table  5-3.  The  results  of  this  calculation  showed  a  small 
negative  gain  of  about  Oq=-.0025  cm  ^  throughout  the  entire  nozzle  for  the 
CC>2(00o  1)  -♦  CO2(10°0)  transition. 

In  conclusion,  the  analysis  of  the  SPGG  device  showed  that  any 
substantial  departure  from  chemical  equilibrium  of  the  flow  in  the  combustion 
and  mixing  zone  would  not  equilibrate  by  the  time  the  flow  had  reached  the 
nozzle  exit  plane.  The  analysis  also  showed  that  there  was  not  a  population 
inversion  hetween  the  (00°  1)  and  (10°0)  vibrational  levels  of  C02  despite 
the  large  departure  from  chemical  equilibrium.  It  is  felt  that  the  intermediate 
levels  of  pressure  and  temperature  (lower  than  most  rocket  motors  and  higher 
than  known  laser  devices)  accounts  for  this  effect.  In  other  words  the  pro¬ 
duction  rates  at  these  pressures  and  temperature  are  not  fast  enough  to  produce 
equilibrium  but  are  too  fast  to  allow  an  inversion. 
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TABLE  5-3 


EXCITED  STATE  C02  REACTION  SET 


RATES* 

A 

N 

B 

(KCALS) 

. 31xl0~4 

-4.9 

0 

.215 

-4.2 

0 

.414 

-4.18 

0 

4.1 

-3.6 

0 

.039 

-4.4 

0 

1.08 

-4.0 

0 

24.4 

-3.28 

0 

2.25 

-3.84 

0 

2.25 

-3.84 

0 

1.365xl0+11 

-  .5 

0 

2.73xl012 

-  .5 

0 

6.64 

-2.59 

0 

.0708 

-3.05 

0 

12 

.43x10  * 

-  .5 

0 

.43xl010 

-  .5 

0 

.43xl012 

-  .5 

0 

12 

.43x10^ 

-  .5 

0 

12 

.43x10** 

-  .5 

0 

C02(110)+IVUC02(020)+M** 

co2(iio)+m-*co2(ioo)+m 

C02  (030)+M-*C02  (020)+M 
C02(030)+M-*C02(100)+M 
C02  (100)+M-»C02(010)+M 


C02(020)+M-*C02(010)+M 

C02(010)+M-*C02(000)+M 

C02(001)+M-*C02(110)+M 

C02(001)+M-»C02(030)+M 

CO9(I10)+M-»COo(030)+M 

CQ2(100)+M-»C02(020)+M 


C02(001)+C02(000) 


►C02(100)+C02(010) 

>C02(020)+C02(010) 


C02(001)+C02(000) 
C02(110)+C02(000) 
C02(030)+C02(000) 
C02(030)+C02(000)-» 
C02(100)+C02(000)-*2  CO2(010) 


C02(100)+C02(010) 

C02(100)+C02(010) 

C02(020)+C02(010) 


C02(020)+C02(000). 


►  2  CO2(010) 


*  Rates  are  for  left  to  right  and  are  of  the  form  AT  N  EXP(-B/RT) 
**A1I  third  body  reaction  rate  ratios  were  taken  as  1.0 
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APPENDIX  A  -  MIXING  MODELS 


This  appendix  describes  two  mixing  models  used  under  this  contract. 
Section  A-l  describes  a  stirred  reactor  model  and  Section  A~2  describes  a 
model  for  the  mixing  in  a  shear  layer  between  two  co-flowing  streams. 

A-l  Stirred  Reactor  Model 

In  a  stirred  reactor  model  reactants  enter  a  reactor  and  are  assumed 
to  be  instantly  and  completely  mixed  with  the  products  in  the  reactor.  The 
combustion  products  extracted  from  the  reactor  are  thus  the  same  as  those 
in  the  reactor  and  only  depend  on  the  input  reactants  and  the  average  stay 
time  in  the  reactor.  This  moc.el  is  particularly  useful  since  it  represents 
the  limiting  case  where  fluid  physics  processes  are  so  fast  as  to  make  the 
chemistry  clearly  the  rate  controlling  step 

In  the  steady  state  process  the  amount  of  mass  in  the  reactor  does 
not  change  and  hence  continuity  becomes 


(A-l.l) 


The  mass  flow  into  the  reactor  is  known  and  consists  of  species 
whose  mass  fractions  are 


ck  in  *  Dk ^  D'in 

The  input  enthalpy/unit  mass ,  h.  ,  is  known,  i.e. 

^2 

hin  =  2(ckhk>in+1 r 

and  the  stirred  reactor  will  be  operating  at  a  known  pressure,  P.  The 
flow  into  the  reactor  will  be  the  sum  of  the  flow  rate  for  each  species 
component 


£ 


(A- 1.3) 
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eactor  which  depends  on  the  How  ie,e  tkro^h  the  reactor  end  the  anenn-4 

ness  in  the  reactor.  The  change  in  species  concetdration  doe  to  ch«  cel 

reactions  -hich  occur  -ithin  a  tine,  r .  is  knoro  to  he  (see  *^7 
equation  (1)  of  Ref.  A-l.)  ^  2.0. 
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Since  the  mixture  is  instantly  stored.  „i„  be  constant  throughout  the 
ctor.  See  reference  A-l  for  a  core  complete  discussion  of  the  checica- 
rate  equations.  Specie  continuity  thus  retires  .ha,  for  each  specie! 
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of  the  reactor  such  that  the  above  system  of  non-linear  sitrul- 
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SHEAR  LAYER 


In  actuality  the  phenomena  of  mixing  in  the  shear  layer  is  highly  two 
dimensional  and  any  simple  one  dimensional  model  should  be  used  with 
care,  e.g.,  see  Ref.  A- 2. 

References  A-3  and  A-4  have  shown  that  the  growth  of  the  shear 
layer  between  2  planar  jets,  co-axial  jets,  and  co-flowing  streams  is 
governed  by  the  growth  law 

b  =  kxN  (A- 2.1) 

where 

b  =  width  of  the  shear  layer 

k  -  constant 

x  =  distance 

N  =  1/2  for  laminar  flow 

=  1  for  turbulent  flow 
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mm 


n— ri  ii  1 1  —i 


In  order  to  approximate  shear  layer  mixing  we  consider  3  stream 
tubes,  the  two  primary  streams  and  one  stream  tube  which  contains  a  small 
fraction  of  the  total  flow  which  is  considered  initially  mixed.  Then  Equation 
(A-2.1)  can  be  rewritten  as 


b  =  b0+k(x-xQ)N  (A-2.2) 

where 

bg  =  width  of  the  initially  mixed  stream 
tube,  i.e.  ,  the  width  of  the  shear 
layer  at  Xg. 

If  it  is  assumed  that  mass  is  entrained  into  the  shear  layer  directly 
proportional  to  the  ratio  of  the  area  of  the  shear  layer  to  the  area  of  the 
primary  streams,  then  the  problem  of  mass,  energy,  momentum,  and 
species  addition  to  the  shear  layer  reduces  to  the  simple  geometric  relation¬ 
ships  given  below: 


planar  flow 
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g(l-fl) 
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kN  (x-Xg) 


N-l 


(r4  -  b.  J 


iO 
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where  r.  =  height  of  the  ith  primary 
stream 

b.g  ~  height  of  the  shear  layer 

protruding  into  the  ith  primary 
stream  at  Xg 

a  -  fraction  of  the  total  mass 

flux  in  the  ith  primary  stream 


fraction  of  the  total  mass 
flux  assigned  initially  to  the 
shear  layer 


co-axial  case  (inner  stream  tube) 


m 


g(i-8)  2(ri-bi)  kN(*-xo’ 


N-l 


b10>* 
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where  r.  =  radius  of  the  inner 
stream  tube 

=  height  of  the  shear  layer 
proturding  into  the  inner 
stream 


co-axial  case  (outer  stream  tube) 


m 


a  (1-fl) 
8 


2(r1  +  b2)  kN  (x  -  xQ)N  1 

r'>  —  (r  +  b 

2  vrl  °20 


(A-2 .5) 


where  r0  =  radius  of  the  outer 

^  stream  tube 

b2  =  height  of  the  shear  layer 
protruding  into  the  outer 
stream 


for  all  cases 

M  =  mv  (A-2. 6) 

H  =  ~  (h  +  y  )  (A-2. 7) 

S.  =  ftic.  (A-2. 8) 

li 

In  the  above  equations  m,  M,  H,  and  are  the  contribution  to  the  mass, 
momentum,  energy,  and  species  addition  rates  respectively  of  the  shear 
layer  from  one  of  the  primary  streams.  For  usage  of  this  option  in  the  computer 
program,  the  reader  should  refer  to  Ref.  A-l. 
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APPENDIX  B 
SUMMARY  OF  RESULTS 

FOR  H/F  and  D/F  IGNITION  COMPUTATIONS 


In  order  to  demonstrate  the  use  of  the  GKAP  computer  code  as  applied 
to  a  problem  of  current  interest  to  AFRPL  a  study  of  the  f^/^/He  ignition  was 
performed.  The  study  included 

Run  #1  :  A  basic  H2/F2/He  (ground  state  reaction  mechanism) 
ignition  calculation 

Run  #2  :  The  effect  of  Deuterum  substitution  on  the  ignition 

Run  #2  :  The  effect  of  considering  excited  state  HF  (vibrational 
level  s  5)  during  the  ignition  calculation. 

Results  of  these  computations  are  summarized  below:  * 

Run  $  1  :  Adiabatic  ignition  of  f^/l^/He  at  constant  pressure.  The 

initial  concentrations  and  conditions  as  well  as  the  reaction 
mechanism  are  presented  in  Table  1.  The  calculation  was 
continued  to  "Equilibrium"  (at  equilibrium  derivatives  of 
T  &  C.  with  respect  to  time  will  be  zero) . 

Table  2  presents  a  comparison  of  an  ideal  Pressure-Enthalpy 
Equilibrium  calculation  with  the  results  of  the  kinetic  cal¬ 
culation  after  ~  7.91  millisec.  As  can  be  observed  from 
Table  2,  the  kinetics  calculation  predicts  the  equilibrium 
conditions  to  within  0.1%. 

Figure  1  presents  the  temperature  profiles  vs.  time  for  all 
three  calculations.  As  can  be  observed  from  Figure  1,  the 
ground  state  H2/F2/He  ignition  model  for  initial  T=300CK 
predicts  a  4.625  millisec  ignition  delay. 

An  interesting  aspect  of  all  three  calculations  may  also  be 
observed  in  Figure  1.  Although  the  final  "equilibrium"  Temperature 
is  s-  1760°K,  during  17  p  sec.  the  temperature  exceeds  2000°  K 
and  has  a  maximum  of  ^  3000° K.  This  temperature  overshoot  is 
due  to  a  lag  in  the  dissociation  of  the  excess  F2.  Figure  2 
presents  the  reaction  net  production  rates  during  the  actual 
ignition.  As  can  be  observed  from  Figure  2,  the  F2  dissociation 


♦Approximately  one  minute  of  CDC  6600  computer  time  is  required  for  each  run. 


lags  the  actual  ignition  and  requires  5-10  p  sec.  to  absorb  the 
excess  energy  via  dissociation.  A  sample  calculation  showing 
the  mechanism  and  the  results  of  simple  hand  calculations 
further  illustrates  the  actual  temperature  overshoot  mechanism 
and  is 

Run  #2  :  Adiabatic  ignition,  at  constant  pressure  for  the  same  conditions 

as  Run  #1  with  Deuterium  exchanged  for  the  Hydrogen  (see  Table  3). 

Figure  1  presents  the  temperature  profile  vs.  time  for  this 
calculation.  The  observed  ignition  delay  is  32.798  millisec. 

The  computer  run  max  paged  after  32.81  millisec  and  was  not 
restarted.  This  ignition  delay  may  be  compared  with  4.625 
millisec  for  the  corresponding  H2  condition,  i.e.  ignition 
delay  ratio: 

TD2  _  32.793  _  _  . 

Th2  4.625 

It  is  interesting  to  note  that  the  reaction  mechanism  for  the  F2 
dissociation  remains  unchanged  from  Run  #1  but  the  "ignition" 
reactions  have  been  slowed  via  D  substitution  for  H.  Therefore 
the  relative  rate  at  which  F2  may  dissociate  and  absorb  excess 
energy  due  to  ignition,  should  increase.  This  should  then  re¬ 
duce  the  temperature  overshoot.  In  fact  a  reduction  of  227° K 
in  peak  temperature  below  the  Run  #1  peak  temperature  is  calcula¬ 
ted.  i 

Run  #3  :  Adiabatic,  constant  pressure  ignition  for  ^A^/He  including 

vibrationally  excited  HF  ( v  s5) .  The  initial  mole  fractions  and 
and  reaction  mechanism  are  presented  in  Table  4. 

Figure  1  presents  the  temperature  profile  vs.  time  for  this  cal¬ 
culation.  The  observed  ignition  delay  is  6.0  millisec  as  com¬ 
pared  with  a  4.625  millisec  delay  for  the  ground  state  mechanism. 


...  J 
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The  observed  difference  in  ignition  delay  may  be  caused  by 
the  HF  exothermic  production  producing  vibrationally  excited 
HF  with  resultant  lesser  energy  channeled  into  thermal  energy. 
This  would  lengthen  ignition  de’~y.  However,  differences  in 
the  reaction  rate  constants  make  this  unclear*  For  example: 

a)  ZH^  H2  Run  #1  k  =  1.0  10"d  T’1 

2H  if  H9  Run  #3  k  =  2.5  1018  T*1 

b)  H  +  F2  if  HF  +  F  Run  #1  k  =  1.21  •  10  e  v  RT  ' 

H  +  F2  ^  HF  +  F  Run  #3  =S  all  reax  for  HF(0-5) 

( _  2.400x 

Run  #3  Z)  -  .8  10 14  e  RT 


No  attempt  was  made  to  determine  effects  of  changes  in 
reaction  rates  on  ignition  delay. 

Figure  3,  presents  the  HF(0-5)  mole  fractions  vs.  time  for 
Run  #3.  As  can  be  observed  from  Figure  3,  an  overproduction 
of  all  vibrational  levels  does  occur,  in  the  same  manner  as 
observed  for  the  ground  state  reaction  mechanism  for  Run  #1. 

One  should  note  that  the  maximum  temperature  observed  in 
Run  #3  =  3006°  K  compared  to  2997°K  for  Run  #1.  This  demonstrates 
that  the  temperature  overshoot  is  nearly  independent  of  the 
vibrational  reaction  mechanism. 


♦Reaction  rates  for  Run  #3  were  taken  from  Kerber  &  Whitier  whereas  for  Run  #1 
and  #2  rates  used  were  supplied  by  AFRPL  for  the  purpose  of  duplicating  results 
obtained  with  the  Aerospace  NEST  Computer  Program. 

Reference:  Kerber,  R.L.  and  Whittier  J.  S. ,  "The  Effects  of  Various  Parameters 
on  the  Behavior  of  the  Pulsed  HF  Chemical  Faser,  Aerospace  Report  No.  TOR- 
0059  (6753-10)-!,  70  December  15. 
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In  addition  to  the  three  computer  runs  discussed  above,  a  series  of 
computer  runs  were  also  made  for  the  purpose  of  examining  the  effect  on 
HzA'z^q  ignition  delay,  r,  with  varying  initial  temperature,  T^.  Results 
of  these  runs,  which  repeat  Run  #1  with 

Tq  =  300,  400,  500,  1000°K 

are  shown  in  Figure  4,  which  gives  Temperature  vs.  time  for  the  above  initial 
temperature  vaiues.  Shown  in  Figure  5  is  a  plot  of  ignition  delay,  r,  vs. 

1000Aq°K  which  can  be  seen  to  be  nearly  linear. 

Simplified  Explanation  of  H2/F2  ^9nition  Results 

In  order  to  check  the  results  of  the  computer  calculations  of  Run  #1  it 
can  be  assumed  that  the  H2/F2  ignition  proceeds  in  essentially  two  steps: 

Step  1) 

H2  +  F  dHF  +  H 
F2  +H  dHF  +  F 

Step  2) 

F2?2F 

The  Step  1  reactions  proceed  much  faster  than  Step  2  which  is  the  dissociation 
of  F  by  third  body  collisions.  Assuming  the  above  two  step  process,  the  temperature 
overshoot  can  be  explai.  'd  by  the  following  calculation: 

For  Step  1  the  result  can  be  written  as 

M(H2)  +  M(F2)  +  N(F2)  d  2M(HF)  +  N(F2) 

For  Step  2  the  equilibrium  result  ran  be  written  as 

N(F2)  d  (N-6)  (F2)  +  2  6(F) 
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Discussion  of  Step  1 


Initial  conditions  were  specified  in  mole  fractions  of  and  ?2  such  that 

M  =  .l59  M  =  .159  mole  H2 

N  =  .302  M  +N  =  .461  mole  F2 

Thus  Step  1  gives  a  mole  fraction  of  HF  of 

HF  =  .316 

i.e. 

.159  H2  +  .159  F2  +  .302  F2  =  .318  HF  +  .302  F2 

The  average  values  for  and  Mw  of  the  mixture  over  Step  1  are 

C  =  .36  cal/g-°  K 
p 

M  =19  g/mole 
w  3 

so  that  since  1  mole  of  HF  has  a  heat  of  formation  of 


AH°  =  -  64800  cal/mole 

then  the  temperature  rise  of  the  mixture  after  Step  1  (producing  .318  moles  of  HF) 
would  be  approximately 

aT  = 

where 

AH  =  .318  *  64800  cal/mole 
0^  =  .36  *  19  cal/mole  0  K 

so  that 

T  _  20600  cal/mole 
~  6  T8  4~cal/mole  u 

AT  =  3010°  K 

Run  #1  gives  a  temperature  rise  of  2700°  K  after  4.625  millisec.  During  this  time 

2  3 

the  species  production  rates  of  Step  1  are  approximately  10  to  10  faster  than 
the  rate  for  Step  2. 
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Discussion  of  Step  2 


The  kinetic  calculations  of  Run  #1  proceed  to  chemical  equilibrium  and  give 
6=  .228 

so  that  Step  2  gives 

.302  (F2)  =  . 074(F2)  +  .456(F) 

The  average  values  for  C  and  M  of  the  mixture  over  Step  2  are 

p  w 

C  =.385 
P 

Mw  =  16.8 

so  that  since  1  mole  of  F  has  a  heat  of  formation  of 
AHj  =  16500  cal/mole 

then  the  temperature  change  of  the  mixture  after  Step  2  (producing  .456  moles  of  F) 
would  be  approximately 

=  AiL 


where 


so  that 


AT  =J± 
C 


ah  = -.456  *  16500  cal/mole 
Cp  =  .  385  *  16.8  cal/mole°K 


T  _  -7520  cal/mole 
A  6.47  cal/mole"  K 


AT  =  -  1160°  K 

Run  #1  gives  a  temperature  drop  during  Step  2  of  1240°  K. 

The  above  simplified  analysis  shows  that  the  kinetic  calculations  of 
Run  #1  give  results  which  are  qualitatively  correct,  including  the  observed  temp¬ 
erature  overshoot  during  ignition. 


* 
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TABLE  1 

H2/T2/HE  Ignition  Calculation  Constant  Pressure 

Initial  Conditions 

P  =  0.42877  ATM 

T  =  300°  K 

Species  Mole  Fractions 


H 

0.0 

HE 

0.38 

F 

2.0E-7 

HF 

0.0 

H2 

0.159 

F2 

0.461 

Reactions  (Rate  Constant  supplied  by  AFRPL)* 


H  +  H 

F2 

-* 

-♦ 

H2 

F  +  F 

n  n 

rH  C 

.V  ^ 

1-0,  11,18  T'‘  35.300, 

5.0  •  1013  e  RT 

HF 

-♦ 

H  +  F 

£ 

k3  = 

.  ,  134.100, 

4.7  •  10  0  T  e  V  RT  ' 

F2  +  H 

-♦ 

HF  +  F 

K  = 

14  _(  2-400, 

1.21  •  10  e  1  RT  ' 

H2  +  F 

-» 

HF  +  H 

4 

k5  ‘ 

14  _  (  1*  600, 

1.60  •  10  e  '  RT  ' 

Reactions  1,2,3  contain  a  generalized  third  body. 


♦units  (cc,  mole,  sec,  °K,  kcal) 
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TABLE  2 
RUN  #1 

Comparison  of  Equilibrium  Calculation  and  Kinetic 
Calculation  at  Time  =7.91  millisec.* 


Equilibrium 

Kinetic 

Difference 


Mole  Fractions 

T°K 

H 

HE 

F 

HF 

H2 

F2 

1607 

0 

.2932 

.4570 

.2453 

0 

.0045 

1603 

10-13 

.2933 

.4562 

.2455 

i<r1B 

.0050 

4 

0 

-  .0001 

+  .0008 

-.0002 

0 

-.0005 

*The  equilibrium  results  were  computed  using  the  NASA/Lewis  Research 
Center  chemical  equilibrium  computer  program.  The  results  for  Run  #1 
shown  here  differ  from  1  gure  1  in  that  the  heat  of  formation  of  F  was  taken 
as  16.5  rather  than  18. 8K  cal/mole  so  as  to  be  consistent  with  the  value  used 
in  the  NASA/LRC  program. 

The  equilibrium  temperature  for  the  kinetic  calculation  has  been 
corrected  from  1598°K  to  1603°  K  by  reducing  the  velocity  from  19. 1  ft/sec 
to  zero  as  follows: 

19 , 1  ft  =>1.6  cal 

sec  - a—  4.9°  K. 

•  33  cal 
g  -°k 

Kin  Temp,  corrected  to  0  velocity  =  1598.0 

+  4.9 
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TABLE  3 
RUN  v2 


D2/F2/HE  Ignition  Calculation  Constant  Pressure 

Initial  Conditions 

P  =  0.42877  ATM 

T  =  300°  K 

Species  Mole  Fractions 


D 

0.0 

HE 

0.38 

F 

2.0E-7 

DF 

0.0 

D2 

0.159 

F2 

0.461 

Reactions  (Rate  Constant  supplied  by  AFRPL)* 


1. 

D  +  D 

-4 

D2 

kj  =  0.3  •  1018  T 

-1 

2. 

3. 

Same  as  Run  #1 

DF  -*  D  +  F 

k»  =  0.47  1019  T“ 

, 136.800  . 

1  V  RT  1 

e 

4. 

F2  +  D 

-4 

DF  +  F 

O 

k.  =0.2  1014  e 

, 2.400, 
t  RT  ’ 

5. 

D2  +  F 

-4 

DF  +  H 

k5  =0.16  1014  e 

,  1.600. 
i  RT  ' 

reactions  1,  2,  3  contain  a  generalized  third  body 


♦units  (cc,  mole,  sec,  °K,  kcal) 
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■&mz  -t 


Sites  ere  resist  so  left  fran  £<e±*r  S  Trritaer 


REACTIONS 
*•2  *  2*h 

H2/F2  ttClTEfi 

#53  wF  T3  5 
•V1*!.®*  3*5. 

1*2*3 

2*F  =  ?2 

!4*5*5*  -8*27.3* 

4> 

«*F  s  l-F(*> 

.4*2.5*  3*23*.’i* 

s 

H*F  *  *“F(1) 

i»*5.5*22i 

4*2.5.  3*122-*"% 

* 

*F(o;  *  hf if i 

A=1.51E7* 

4*-'i.*5S*  ■**!>•  I* 

i*.3$ 

rtFiOJ  =  fFin 

A=l.r*£-2. 

iw*-3.t5«  ®*5.5* 

it 

**F(H  =  »*F  121 

i*2.=^Ti 

4= -I.455.  5*5.5* 

17.2 

**F  (1 1  =  HF  (2> 

A=2.3*E-2. 

V=-3.65,  5*5.5* 

1? 

hF(21  *  Hr  (31 

A*4.S*E7, 

4,*— 1 .455  *3*®.®* 

25*21 

*«F  (21  *  r*F  (31 

a=3.57e-2. 

V»*J.‘S*  3*5.5* 

22 

h:F(3J  *  hF(*1 

A=6- 1*F7 . 

flw-1  •  455 ,3*5.®* 

23*24 

hF (31  *  HFI41 

X=4. ?»£-?. 

4?— 3. *5*  5*5.5* 

23 

hF(*1  =  HF151 

IsT.SPt?. 

.4=- 1.455*^5*5.®* 

2**27 

H?{41  *  HF 15 > 

4*5.  2. 

fc*-3*fe5*  8tef.®* 

25 

h *F  -  hF(21 

A=5.ffit-Z2* 

14=2  «  ®  •  =•*  111-  5*4 

.*7 

H*F  =  HF  (31 

A=3.Ct22* 

4*2.4,  3=  m 

E^’D  THR  REAX 
hF  (01  *  H  »  F 

*  H2* 

J=7.©-12» 

ft,1*©.®*  3*1.71* 

7 

hF  ( i  1  ♦  h  3  F 

*  H2* 

A=J4. J£12. 

4*©.C*  3*1.71* 

$ 

MF(?»  ♦  h  3  F 

*  1*2  • 

*=7.e=-i3. 

4=5.5,  3*1 .71* 

9 

F  *  i*2  3  HF  (31 

*H 

f 

A*3. 0-1?, 

\«-i5i  3*5.0* 

15 

F  ♦  H2  =  HF (41 

*H 

• 

**3-0^12* 

3=®.S* 

11 

F  ♦  H2  =  HF (51 

*H 

f 

A33.CL12. 

3=0.0* 

12 

HF(01*FFt2l32*HF(ll 

,A32.  =.4, 

4=-2.2*  3*©.o* 

32 

hF(11*hF13)*2*hFC21 *A=2.ot4. 

I4=— 2  .2  *  3*5 .©• 

29 

HF(2)*HFt4)=2*HF(3l  .A=2.0t4, 


\*-c •  4Q  , 


hF(3)*7FIS1*2*HF(4) «A*2. 

Ct4,  3 S=- 

2.4Q,  tf— 0 

•  Si  41 

hF  I0>**-F(31*HF(1  l*HF(2i  * 

A=6.?5El* 

N=— 2.5* 

0=  ' .  0  * 

45 

HF  n  1  *t-F  141  =HF  (2)  ♦hF  (3l  * 

A=6.25E1* 

M*-2 ,o» 

0*6.©# 

46 

HF  (2l4i-Fl5)aHFl3l*HF(4)  . 

A*6.?5E1* 

K=— 2  a  8  • 

5*0.0* 

47 

hF(0) 4RF ( 4 1 sHF (11 *HF (31 « 

A=3.oE-3* 

N=-3.R» 

3*0.®* 

51 

hF(11*I-F(5)*HFI21*HF(4)  « 

A=3.0E-3* 

0*v.3* 

52 

hF(0)  ♦  F  *  h  ♦  F2, 

A=5.0E12* 

N=0.0* 

13*2 .4©  , 

65 

HF(1)  ♦  F  *  H  ♦  F2i 

A=5.0E12* 

«=0.0* 

0*2. 4Q, 

51 

hF(2)  *  F  *  h  ♦  F2* 

A=in.oE12* 

N*0.0* 

3*2.40* 

62 

HFC31  ♦  F  *  H  ♦  F2» 

A*lo.oE12* 

N*0.0* 

0*2 . 40  * 

63 

hF ( 4 1  ♦  F  a  H  ♦  F2 • 

A=2o.oE12, 

N*0.0* 

0*2.40* 

64 

hF(5)  ♦  F  *  h  *  F2* 

A*3o.oE12* 

N* 0*0* 

5*2.40* 

65 

LAST  hEax 


third  bcoy  rate  ratios 

SPECIES  h2*  4*1 • 0 • 1 « 0 « 1 • 0* 1 » 0, ] , 0 *  1 » 1 * 1  * 

SPECIES  Hi  8*3*1*  0  »  1  *0  •  1  •  0 »  1 , 0*  1  * 0  *  1  *  1  *  1  * 

SoECIES  E*  .4*3*1*  0 *  1  * 0 • 1 • 0 *  1  * 0* 1 » 0 *  1 • 1  *  1  * 
SPECIES  f2  *  >4 *  1 « 1 • 1  *  1  * 0 • 1 « o, 1 « 0* 1  * 0 1 1  * 0 *  1  *  1  * 

SPECIES  HF(0) *.4, 1*1, 1*1* 0*1, 0,1* 0*1, 0*1, 0*1*1, 
SPECIES  HF(l)«.4*ltl«l*l*0*l(0,l*0*ltP>lf0*l*l* 
SPECIES  HF (2) *«4*l*l*]fl*0*l»0*lt0*l*0*l*0tl«i* 
SPECIES  hF (3) *,4*l*l*l*l*0*l*0«lt0*l*0»lt0*ltl« 
SPECIES  HF (4),. 4, 1,1, 1,1, 0*1, 0,1*0* 1*0*1, 0,1*1* 
SPECIES  HF (5) ,.4*1* 1*1*1* 0*1*0. 1*0»1*0*1*0*1*1. 
SPECIES  HE* ,4*3*1 »  0 *  1 • 0 *  1 • 0 *  1 , o » 1  * 0 » 1  * 1  *  1  * 

LAST  CARD 


species  *cle  fr.sctio’v'S 


H2 

H 

F 

F2 

HF  (01 
HF(11 
HE  (2  1 
HF  (31 
HF  (41 
HF  (51 
HE 


.155 

0.0 

o.onso  02 

.*51 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

.380 
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See  Reaction  Table  1 


See  Reaction  Table  2 


H2A2 
ground 
state 
Run  #1 


VK 

1000 

T 

T 

300 

3.33 

0 

400 

2.5 

1.71 

500 

2. 

2.88 

1000 

1. 

4.62 

5 


Ignition  Delay 

4 

T 

m  Hi  sec. 


3 


2 


1 


0 


I 


3  POO 

V* 


'i-jrre  -5;  Tbe  ctf  TesrpErat-ar^cc  a^nilirr-  Delay 

fcr  H2/T2'Be  -  Ts!0  =  2*  W“'T 


?-  2  -5 


0 


1.0 


2.0 


3.0 


.0 


TIME  (MILLI5EC.) 

Ti—jje  S:  HZ/T2/HL  Ignition  I  =  300.  490, 
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